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We discuss the phase structure of chiral Yukawa models in the mean-held 
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phase transition lines near a critical point. Our results contrast with some 
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I. INTRODUCTION 



Chiral Yukawa models on the lattice are interacting fermion-Higgs field theories with 
a chiral symmetry, which are useful for studying the Standard Model on the lattice in 
the absence of gauge couplings [|l|. The phase structure of such models, in particular the 
location of fixed points and second order phase transition lines, is of interest for the definition 
of corresponding continuum field theories. 

The action of a chiral Yukawa model is of the generic form 

+ yY,*x(PL®+ + PR$ x )y x (1) 

x 

+ fermion kinetic term + Higgs potential, 

the detailed definition of these terms varying from one model to the other. With 'triviality' 
in mind, we take the bare Higgs self-coupling equal to oo, which fixes the radial mode 
of the Higgs field $. The phase diagram also depends on the details of the model, but 
the phase structure in the small to intermediate-?/ region is similar in most models. We 
shall concentrate on this part of the phase diagram (see Fig. H), which is sufficient for our 
purposes, but the considerations presented here are also valid for critical points at large-?/, 
where applicable. There are paramagnetic (PM), ferromagnetic (FM), antiferromagnetic 
(AM) and possibly ferrimagnetic (FI) phases, characterized by the vacuum expectation 
values of the spatial average and staggered average (this means that fields at odd sites 
contribute with a minus sign) of the Higgs field and the fermion condensate. In the context 
of chiral Yukawa models, phases where both the Higgs field and its staggered analogue have 
non-zero expectation value are known as FI phases 0. 

In recent years, various kinds of Yukawa models, based on different lattice fermion actions 
and with different symmetry groups, have been studied both using mean-field methods and 
numerically . Here we shall concentrate on the mean-field approach || . Because of the 
difficulty of dealing with the fermion determinant exactly, mean-field techniques are often 
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applied in combination with small or large-?/ expansions. In Ref. |4| this type of analysis 
was extended to FI phases. Although FI phases start off at values of y where one might 
question the validity of the small-?/ expansion involved, it was argued that the real expansion 
parameter is the quantity y(^f^f), which remains small as long as one keeps close to the PM 
phase. Furthermore, the phase structure of the models studied, obtained from numerical 
simulation, appeared to be well-described by the mean-field results. 

In a recent publication ||, mean- field calculations were presented for a general class of 
chiral Yukawa models with different symmetry groups and lattice fermion actions. (The 
details of this approach differ somewhat from other applications of mean-field methods in 
the truncation of the fermion determinant, but this is of no importance here.) The authors 
claim that, in general, if the PM-FM and PM-AM phase transition lines intersect at a 
point A (cf. Fig. [I]), then they must necessarily continue beyond the point A with smooth 
first derivatives at A. There would always be a FI phase, with the PM-FM line continuing 
as an AM-FI phase transition line and the PM-AM line as an FM-FI line. (Everything 
within the mean-field approximation.) This is in contrast with the results of Refs. [f|||, 
where discontinuities were observed in the slopes of the phase transition lines at A, with 
indications in one model || that a FI phase is absent at all. (Again, within the mean-field 
approximation. ) 

These discrepancies motivated us to carry out the present study, which may however have 
a wider applicability. We analyse the phase structure around the point A from a general 
point of view, in the mean-field approximation. It is found that the first derivatives of the 
second order phase transition lines, assuming that they intersect, are in general discontinuous 
at A, as in Fig. [I], and no conclusion can be drawn a priori about the existence of a FI phase. 
After presenting a general demonstration in Sec. 0, we illustrate the results with a simple 
example in Sec. [TIT]. For definiteness we keep in mind a phase structure as in Fig. |I[ 
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II. PHASE STRUCTURE AT A 



In the mean-field approximation, the phase of the system at a point (y, k) is determined 
by minimizing the free energy F with respect to a number of mean fields h l and staggered 
mean fields h\ (i — 1, . . . ,N), collectively denoted as h, h s . (F is a function of h and h s 
with coefficients depending on y and k.) This gives the mean- field equations 

OF OF 

IK = °- as; = < 2 > 

A solution of these equations corresponds to a (local) minimum of F if the matrix F" of 
second derivatives at the solution is positive definite. A second order phase transition occurs 
when one of the eigenvalues of this matrix goes through zero: a negative mode develops, 
destabilizing the original solution and replacing it by one with a lower free energy, belonging 
to a different phase. Therefore, the condition for a second order phase transition line is 
given by the additional equation 

detF" = 0. (3) 
If there is only one mean field h and one staggered mean field h s this becomes^ 

0. (4) 



d 2 F d 2 F ( d 2 F \ 2 



dh 2 dh 2 \ dh dh s J 

In the chiral Yukawa models under consideration, the free energy F is symmetric under 
a sign change of all the fields h l at the same time or all the fields h\ at the same time. This 
implies the absence of terms containing odd powers of h or h s in the Taylor expansion of F, 
such as hh s . Furthermore, at the four phase transition lines under study here, either all the 
h % = or all the h\ = or both. Hence d 2 F/dh l dh{ = there and Eq. (J3|) simplifies, 



1 There is a square missing in the corresponding formula (17) of Ref. 
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Which of the two factors in this expression is zero depends on which phase transition line 
one considers; at the intersection point A both factors are zero. 

For definiteness we shall focus on the PM-FM and AM-FI transitions. Then condition 
(H) becomes 

det|^(/i = 0A) = 0, (6) 

since in this case the unstable mode is in the direction of the h fields. Along the PM-FM 
line, h s — and the solution to Eq. ([]) is readily obtained. Along the AM-FI line, however, 
h s ^ 0. One has to solve Eqs. @) and @ simultaneously, to find h s (y, k) and the function 
k c (y) parametrizing the line. This procedure was applied numerically for phase diagrams of 
chiral Yukawa models based on the Zaragoza proposal |7j with the 'most local' and 'Roma 
I' || fermion actions f|||. 

Here, we are especially interested in the behaviour of the free energy F in the vicinity of 
the point A. In this region, close to the PM phase, the mean fields are small and one can 
study the phase structure by considering an expansion in h and h s . In Ref. || it is claimed 
that it is sufficient for this purpose to know F up to quadratic terms in h and h s . The 
following calculation of the derivatives of the phase transition lines at A will show, however, 
that quartic terms (or, in their absence, higher order terms) are indispensable; neglecting 
them leads to incorrect statements about a possible FI phase near A. 

The PM-FM and AM-FI phase transition lines are given by a continuous function k c (y) 
whose derivative we wish to determine. Defining 

d 2 F 



f = det 



(7) 

h=0 



dh 2 

we can write Eq. @ for these lines as 

f(yMv)X(vMv))) = Q' (8) 

Taking the derivative of this equation in the direction tangential to the line we find 

o = # = ?l + d J- d ^ ^a^i) (9 ) 

dy dy dk dy dh? s \ dy dk dy J ' 



leading to a slope 



dy \dy dh 2 s dy ) \dk dh 2 s dk )' 1 ' 



On the left hand (PM-FM) side of point A, h 2 s = and Eq. QTIJ) reduces to 

fd_K\ = _V/V. (n) 

\dyJ PM _ FM dy I dk 

On the AM-FI side, however, it is reasonable to assume that h 2 approaches A linearly in 
y—yA and k — kA, in accordance with a mean- field critical exponent of 1/2 for h s in this region 



(cf . the example in Sec. |HI|) . It follows that both the numerator and the denominator of Eq. 
(|10|) receive an additional non-zero contribution on this side, so that dk c /dy is discontinuous 
at A. 

A similar analysis can be carried out for the change of slope between the PM-AM and 
FM-FI lines. We emphasize, however, that the AM-FI and PM-FI lines thus obtained 
(at least, infinitesimally close to A) remain 'candidate phase transition lines' only until 
free energy considerations establish that the phases on both sides of the lines correspond 
to absolute minima of the free energy. The slope of the AM-FI line suggested by these 
calculations may, for example, come out bigger than that of the PM-FI line. (This does in 
fact happen, at least to lowest order in the Yukawa coupling y, in a chiral Yukawa model 
based on the Roma I action ||.) In that case, a comparison of free energy values should 
indicate which of the calculated lines does or do not correspond to a true second order phase 
transition. 

III. A SIMPLE EXAMPLE 

We would like to illustrate these results with a simple model of one mean field h with its 
staggered analogue h s , described by a quartic free energy F, 

F = --ah 2 - -bh 2 + -ch 2 h 2 + — dh A + — eht (12) 
2 2 s 4 s 24 24 s y ' 
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where the parameters a, . . . , e are well-behaved functions of y and k. We assume the presence 
of PM, FM and AM phases as in Fig. [I], i.e., a increases and b decreases with increasing 
k, and a(y, k) and b(y, k) are such that the PM-FM and PM-AM lines meet in a point A. 
Furthermore, for stability of F, we require d > 0, e > 0, c > —\fde/2>. Apart from these 
stability conditions, one can consider this model as the expansion of a general free energy 
in the neighbourhood of the point A up to quartic terms in the fields. 

In those points of the (y, fc)-plane where both a < and b < 0, the mean-field equations 
(§) imply that the system is in a PM phase with h = h s = and free energy normalized to 
zero. If a > or b > we are in one of the broken phases, FM, AM or FI. Condition @ for 
the PM-FM and AM-FI second order phase transition lines becomes 

-a + -ch 2 s =0. (13) 

Along the PM-FM line this becomes simply a = 0, whereas along the AM-FI line we need 
the value of h s which follows from the mean- field equations (0), taken at h = 0, 

-bh s + ^eh 3 s = 0. (14) 

At the point A, both a = and 6 = 0, whereas in general c, d and e are non-zero. Close to 
A we can therefore write, 

a(y,k) = a y (y - y A ) + a k (k - k A ) + h.o., 

b(y,k) = -b y (y - y A ) - b k (k - k A ) + h.o., 

c(y,k) = c + h.o., 

d(y,k) = d + h.o., 

e(y,k) = e + h.o., (15) 

where h.o. stands for higher orders in the y — y A , k — k A expansion. The signs in front of 
afc and bk have been chosen such that both are positive. In terms of these coefficients, the 
slopes of the PM-FM and PM-AM lines at A are (cf. Eq. (|TTJ)) 



7 



i?PM-FM = (^(A)) = (16) 

\ uy J PM-FM ak 

-RpM-AM = _ T~- (17) 



In order to find the AM-FI slope at A we solve Eqs. (13) and (14) to lowest order in y — y A 
and k — k A - We find 

h\ = --{b y {y - y A ) + h(k - k A )) + h.o., (18) 

corresponding to a critical exponent 1/2 for h s as mentioned earlier, and 

ay + 3(co/e )6y 

-Kam— fi — / — rr l iy J 

a fc + 3(c /e )6 fc 

(cf. Eq. ©). 

It is interesting to express Eq. (|T9|) as a 'weighted average' of the PM-FM and PM-AM 
slopes: 

eodfc-RpM-FM + 3c &fc-RpM-AM /or .N 

-Kam-fi — ^ — r • l^UJ 

e a fc + 3c fc 

Similarly, we find for the FM-FI slope 

3Coafci?PM-FM + cWfc-RpM-AM /01 x 

-KFM-FI - ~ — TT • l^J-J 

3c a fc + d b k 

The discontinuities in the slopes are evident. Only for c = the first derivatives of 
these phase transition lines are continuous (ignoring special cases like c 7^ 0,Cq = 0). This 
corresponds to the trivial case that there is no coupling between h and h s , as mentioned 
earlier. Note that in the limit e — > Eq. (0) tends to -Ram— fi — -Rpm-am, while in the 
limit d — > which can be taken simultaneously Eq. ([H]) leads to -Rfm-fi — -Rpm-fm- In 
this case the AM-FI and FM-FI lines would be 'interchanged' compared with the c = case. 
This obviously indicates that at least one of the calculated 'candidate transition lines' does 
not represent a genuine phase transition, and additional free energy considerations must 
determine the real nature of the transitions, as discussed before. 

It is also instructive to consider the AM-FI and FM-FI slopes as a function of Cq. In the 



limit that Co approaches its 'stability lower bound' — y/d e /3, these lines form a 180 degree 



angle. Upon increasing cq this angle decreases, until it vanishes for Co = For still 

larger values of cq, corresponding to strong coupling between h and h s in the free energy, 
the FI phase disappears (at least in a neighbourhood of the point A; higher order terms in 
Eqs. (^) may give rise to a FI phase a little farther out, as in Fig. |), and instead there is 
a first order transition separating the AM and FM phases. 

IV. CONCLUSION 

For determining the phase structure close to the point A where several phases meet, in a 
mean-field approximation, it is not sufficient to expand the free energy up to terms quadratic 
in the mean fields. Quartic terms are crucial for determining the slopes of the transition lines 
enclosing the FI phase, and in general the slopes are discontinuous at A. Such behaviour 
was observed in mean-field studies of phase diagrams of chiral Yukawa models [fH|, where 
the (candidate) transition lines bordering the FI phase were determined numerically. 

These results explain the discrepancy signalled in Sec. |. The authors of Ref. |5j in- 
correctly assumed that it is sufficient to consider the free energy up to quadratic terms in 
the mean fields. As a consequence, the conclusions drawn there about the presence of FI 
phases and the transition lines separating them from the FM and AM phases are incorrect. 
However, the results for the PM-FM and PM-AM lines are not affected. We emphasize 
that these conclusions do not depend on the approximations made in the treatment of the 
fermion determinant in chiral Yukawa models. The actual location and slopes of the var- 
ious phase transition lines do, however, because different approximations lead to different 
'effective' free energies. Hence, improved approximations of the fermion determinant may 
change the conclusion about the existence of a FI phase, apart from the limitations of the 
mean-field approximation. 

This work was supported by EC contracts ERBCHBICT941067 and CHRX-CT92-0051, 
by DGICYT (Spain) and by Accion Integrada Hispano-francesa HF94-150B. 
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FIGURES 

FIG. 1. Generic phase diagram of a chiral Yukawa model in the small-y region. On the vertical 
axis is the hopping parameter k of the Higgs field, on the horizontal axis is the Yukawa coupling 
y between the Higgs field and the fermions \£. 



FIG. 2. Possible (mean-field) phase structure for Co > \/^0 e o/3, with a first order AM-FM 
transition (dashed line) and a FI phase away from the point A. 
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Fig. 1 
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Fig. 2 
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